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Abstract 

Swine enteric coronaviruses are a group of most significant pathogens causing diarrhea in piglets with similar clinical symptoms and 
pathological changes. To develop a simple, rapid, accurate, and high-throughput detection method for diagnosis and differential 
diagnosis on swine enteric coronaviruses, specific primers and probes were designed based on the highly conserved regions of 
transmissible gastroenteritis vims (TGEV) N, porcine epidemic diarrhea vims (PEDV) M, porcine deltacoronavims (PDCoV) M, 
and porcine enteric alphacoronavims (PEAV) N genes respectively A TaqMan-probe-based multiplex real-time RT-qPCR assay was 
developed and optimized to simultaneously detect these swine enteric coronaviruses. The results showed that the limit of detection can 
reach as low as 10 copies in singular real-time RT-qPCR assays and 100 copies in multiplex real-time RT-qPCR assay, with all 
correlation coefficients (R ) at above 0.99, and the amplification efficiency at between 90 and 120%. This multiplex real-time RT-qPCR 
assay demonstrated high sensitivity, extreme specificity, and excellent repeatability. The multiplex real-time RT-qPCR assay was then 
employed to detect the swine enteric coronavims from 354 field diarrheal samples. The results manifested that TGEV and PDCoV were 
the main pathogens in these samples, accompanied by co-infections. This well-established multiplex real-time RT-qPCR assay 
provided a rapid, efficient, specific, and sensitive tool for detection of swine enteric coronavimses. 
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Introduction 

Porcine enteric viruses, the pathogens of viral diarrhea in pigs, 
caused huge economic losses in the swine industry in recent 
years (Zhang et al. 2016). There were more than ten enteric 
viruses discovered from swine gut in the past, including but not 
limited to transmissible gastroenteritis vims (TGEV), porcine 
epidemic diarrhea vims (PEDV), porcine deltacoronavims 


Electronic supplementary material The online version of this article 
(https://doi.org/10.1007/s00253-019-09835-7) contains supplemental 
material, which is available to authorized users. 




Jinquan Wang 

Wangj inquan 163 @ 163. com 




Guangliang Liu 
LiuGuangliangOl @ caas.cn 


College of Veterinary Medicine, Xinjiang Agricultural University, 
311 EastNoda Road, Urumqi 830052, Xinjiang, China 

State Key Laboratory of Veterinaiy Etiological Biology, Lanzhou 
Veterinary Research Institute, Chinese Academy of Agricultural 
Sciences, 1 XuJiaPing, YanChangBu, ChengGuan District, 
Lanzhou 730046, Gansu, China 


(PDCoV), porcine enteric alphacoronavims (PEAV), porcine 
sapovirus (PSaV), porcine norovirus (PNoV), porcine 
teschenvims (PTV), porcine kobuvims (PKV), Seneca Valley 
vims (SW), porcine rotavirus (PRV), porcine reovims (ReoV), 
porcine bocavims (PBoV), porcine astrovims (PAstV), and etc. 

Within these swine enteric viruses, coronavimses are the 
most devastating pathogens responsible for acute diarrhea, 
vomiting, dehydration, and high mortality in neonatal and 
suckling piglets. According to the genetic and antigenic char¬ 
acters, coronavims was divided into four genera: Alpha-, Beta-, 
Gamma-, and Delta-CoV (Woo et al. 2010). During the past 
years, some alphacoronavimses (PEDV, TGEV, and PEAV) 
and a deltacoronavims (PDCoV) emerged or reemerged in 
the pig farms and resulted in severe diarrhea and mortality in 
the early stage of suckling piglets. 

The TGEV and PEDV were the traditional causal agents re¬ 
sponsible for diarrhea in pigs for the past decades. However, the 
variant strains of PEDV were discovered since 2010 and circu¬ 
lating in pig herds thereafter, showing up to 100% death rate for 
piglets younger than 1-week-old (Lyoo et al. 2017; Sun et al. 
2012). Another swine enteric vims, porcine deltacoronavims 
(PDCoV), was firstly discovered from healthy pig herds by a 
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research team of Hong Kong in 2012 when they did a molecular 
epidemiology study in 3137 mammals and 3519 birds (Woo et al. 
2012). Two years later, this vims was reported to cause severe 
diarrhea and/or vomiting and atrophic enteritis in the USA (Wang 
et al. 2014) and later on in China (Song et al. 2015). 
Subsequently, a novel porcine enteric alphacoronavims PEAV 
was discovered in 2017 from some pig farms located in 
Southern China, causing more than 24,000 piglets’ death 
(Gong et al. 2017). This vims was also named as SeACoV 
(Pan et al. 2017) or SADS-CoV (Zhou et al. 2018b) by the other 
research groups. 

All these four swine enteric coronaviruses, causing 
similar clinical symptoms and pathological changes in 
piglets, circulate in the pig herds and result in huge eco¬ 
nomic losses across the world in recent years. To devel¬ 
op a simple, rapid, accurate, and high-throughput detec¬ 
tion method for diagnosis and differential diagnosis on 
swine enteric coronavimses, a TaqMan-probe-based mul¬ 
tiplex real-time RT-qPCR assay was established to simul¬ 
taneously detect TGEV, PEDV, PEAV, and PDCoV from 
the same reaction vial. To our knowledge, this is the first 
multiplex real-time RT-qPCR method for swine enteric 
coronaviruses detection. 

Materials and methods 

Viruses, dinical samples, primers, and probes 

The TGEV was previously propagated and preserved in our 
laboratory. The PEDV and PDCoV were isolated from clinical 
samples and confirmed by conventional PCR and DNA sequenc¬ 
ing (GENEWIZ, Suzhou, China). The nucleoprotein gene of 
PEAV (GenBank access number: MF370205) was synthesized 
from GENEWIZ biotech company (Suzhou, China). 

Field samples were collected from diarrheal piglets between 
2015 and 2018 from Liaoning, Shandong, Chongqing, Shaanxi, 
Ningxia, and Gansu provinces. All samples were stored at - 
80 °C until use. 

The primers and TaqMan probes for real-time qPCR assay 
were designed by the software Beacon Designer 7 (PREMIER 
Biosoft International, Palo Alto, CA, USA). The detailed in¬ 
formation of primers and probes were listed in Table 1 . 

RNA extraction and reverse transcription 

All clinical samples were resuspended with phosphate- 
buffer saline (PBS), vortexed and centrifuged at 
12,000xg at 4 °C for 10 min. An aliquot of 250 pU 
supernatant was applied for total RNA extraction with 
RNAiso reagent (Takara, Dalian, China) following the 
manufacturer’s instruction. The total RNA in RNase- 
free water was reversely transcribed into cDNA using 


RevertAid First Strand cDNA Synthesis Kit (Thermo 
Scientific, Waltham, USA). The cDNA was subjected 
to real-time qPCR analysis with the multiplex RT- 
qPCR method established in this study. 

Construction of recombinant plasmids for standard 
curves 

The M genes of PEDV and PDCoV, and N genes of TGEV and 
PEAV were constmcted into pET-30a(+) vector. The recom¬ 
binant plasmids were linearized by Smal enzyme digestion 
and recovered by PCR purification kit. The purified recombi¬ 
nant plasmids were quantified by spectrophotometric analysis. 
The copy number of recombination plasmids was calculated 
by using the following formula (Huang et al. 2009): 

(6.02 x 10 23 ) x (Xng/ftL x 1(T 9 ) 

Plasmid copies pL = ----A--- 

plasmid length (bp) x 660 

To establish the standard curves for single coronavims, each 
plasmid was diluted in a tenfold series, from 10 copies/pE to 
10 1 copies/pE. For multiplex standard curves, each of the four 
linearization plasmids was adjusted to 4 x 10 9 copies/pE and 
pooled with equal volume to made 1 x 10 9 copies/]LiL of each 
plasmid. The pooled plasmid was then diluted serially by tenfold 
to establish multiplex standard curves. 

Single and multiplex real-time qPCR condition 

All real-time qPCR reaction systems were set to a volume of 
20 pE. For single qPCR amplifying TGEV, PEDV, and 
PDCoV, 10 pE 2 x TransStart Probe qPCR SuperMix 
(TransGene, Beijing, China), 200 nM primers and probe each, 
1 pE plasmid DNA template, and 6.8 pE nuclease-free water 
were pooled and mixed. For PEAV amplification, 10 pE 2 x 
TransStart Probe qPCR SuperMix (TransGene, Beijing, China), 
500 nM each primer, 100 nM probe, 1 pE plasmid DNA tem¬ 
plate, and 6.8 pE nuclease-free water were pooled and mixed. 
All reactions were amplified on a Bio-Rad CFX96™ Real-time 
System (Bio-Rad, Hercules, CA, USA) at 94 °C for 30s, follow¬ 
ed by 40 cycles of 94 °C for 5 s and 60 °C for 30s. 

For reaction system of multiplex real-time PCR, 
10 pE 2 x TransStart Probe qPCR SuperMix combined 
with all primers, probes, templates, and nuclease-free 
water to a final volume of 20 pE. The concentrations 
of each primer and probe of PEDV, PDCoV, TGEV, and 
PEAV were optimized for better outputs. The amplifying 
cycles of multiplex qPCR were carried out as same as 
singular real-time qPCR. The Cq value higher than 35 
was considered negative. All qPCR results were ana¬ 
lyzed by CFX Manager™ software. 
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Table 1 Primers and probes 

designed for real-time RT-qPCR Virus Primer/ Sequence (5'-3') 

probe 


PEDV 

Forward 

Reverse 

Probe 

GATACTTTGGCCTCTTGTGT 

CACAACCGAATGCTATTGACG 

FAM-TTCAGCATCCTTATGGCTTGCATC-TAMRA 

M 

150 

PDCoV 

Forward 

Reverse 

Probe 

ATTTGGACCGCAGTTGACA 

GCCCAGGATATAAAGGTCAG 

Cy5-TAAGAAGGACGCAGTTTTCATTGTG-BHQ2 

M 

92 

TGEV 

Forward 

Reverse 

Probe 

TGCCATGAACAAACCAAC 

GGCACTTTACCATCGAAT 

HEX-TAGCACCACGACTACC AAGC-BHQ1 

N 

81 

PEAV 

Forward 

Reverse 

Probe 

TCTCGGCTTACTCTAAACCC 

CATCCACCATCTCAACCTC 

TexasRed-AAGACCTAAATGCTGATGCCCCA-BHQ2 

N 

150 


Target Size 

gene (bp) 


Sensitivity, specificity, and repeatability analysis 
of multiplex RT-qPCR assay 

To analyze the sensitivity of established multiplex RT-qPCR, 
the linearization standard plasmids prepared above were dilut¬ 
ed tenfold serially to a final concentration between 1.1 x 

H • i 

10 copies/qL and 1.1 x 10 copies/qL in nuclease-free water. 
The diluted standard plasmids were used as templates for real¬ 
time qPCR amplification. 

To estimate the specificity of this established multiplex RT- 
qPCR, standard DNAs, or cDNAs of major swine viruses, 
including porcine astrovirus (PAstV), porcine kobuvirus 
(PKV), type O foot-and-mouth disease virus (FMDV-O), type 
A foot-and-mouth disease virus (FMDV-A), porcine repro¬ 
ductive and respiratory syndrome virus (PRRSV), classical 
swine fever virus (CSFV), porcine rotavirus (PRV), porcine 
circovirus 2(PCV2), and pseudorabies virus (PrV) were used 
as templates for amplification. The nuclease-free water was 
served as negative template control. 

To evaluate its repeatability, tenfold serially diluted stan- 
dard template between 1.1 x 10 copies/qL to 1.1 x 
10 1 copies/qL were used to test the coefficients of variation 
of real-time PCR. For intra-assay repeatability, all samples 
were triplicated. For inter-assay repeatability, the assays were 
repeated three times individually at different locations. 

Clinical sample detection 

A total of 354 fecal samples were collected from diarrheal pig 
farms located in Liaoning, Shandong, Chongqing, Shaanxi, 
Ningxia, and Gansu provinces of China between 2015 and 
2018. All samples were diluted fivefold with sterile 
phosphate-buffered saline (PBS), vortexed, and centrifuged 
at 1847xg at 4 °C for 20 min. The supernatant was collected 
and used to extract viral RNA with Trizol reagent (Invitrogen, 
Carlsbad, CA, USA). The cDNAs were generated by Reverse 
Transcript System (Promega, Madison, WI, USA) using 


extracted total RNA as templates and hexamer random 
primers. All cDNA from clinical samples were measured by 
the multiplex RT-qPCR assay developed in this study. 

Results 

Single real-time RT-qPCR assay for individual virus 

To develop a multiplex real-time RT-qPCR, the single real-time 
RT-qPCRs for detection of the individual virus were firstly 
established with different fluorescence-labeled target probes, for 
details, FAM for PEDV M gene, Cy5 for PDCoV M gene, HEX 
for TGEV N gene, and TexasRed for PEAV N gene. The stan- 
dard curves for each vims were generated using 1.1x10 copies 
to 1.1 x 10 1 copies of tenfold serially diluted linearized plasmids 
conceiving target genes. The results demonstrated that the single 
real-time RT-qPCR assays for each vims were successfully 
established at the limit of detection at approximately 10 copies 
(Fig. 1). All the standard curves showed an excellent correlation 
coefficient and amplification efficacy, for instance, PEDV (R 2 = 
1; Eff% = 95.7), PDCoV (R 2 = 1; Eff% = 100.9), TGEV (R 2 = 
0.999; Eff% = 106.0), and PEAV (R 2 = 0.995; Eff% = 98.9), in¬ 
dicating that the single real-time RT-qPCR for each vims was 
valid and reliable. 

Establishment of multiplex real-time RT-qPCR assay 

To establish the multiplex real-time RT-qPCR method, all primer 
sets, probes, and serially diluted standard plasmids for detection 
of PEDV, PDCoV, TGEV, and PEAV were mixed with 2 x 
TransStart Probe qPCR SuperMix and nuclease-free water. The 
concentrations of each primer and probe were optimized for the 
optimum output. The optimal final concentrations of primers and 
probes were as follows: 300 nM primer and 100 nM probe for 
PEDV, 200 nM primer and 200 nM probe for PDCoV, 500 nM 
primer and 100 nM probe for TGEV, and 500 nM primer and 
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Amplification 




O Standard 
X Unknown 

FAM E=95.7% R A 2=1.000 Slope=3.430 y-int=39.933 


C Amplification 



Standard Curve 



O Standard 
X Unknown 

- HEX E=105.0% R A 2=0.999 Slope=3.187 y-int=36.239 


Fig-1 Establishment of single real-time RT-qPCR assay for the individ¬ 
ual virus, a The amplification curves (top) and a standard curve (bottom) 
for detection of PEDV M gene were generated. The probe was labeled 
with FAM at 5'-end and TAMRA at 3'-end. b The amplification curves 
(top) and a standard curve (bottom) for detection of PDCoV M gene were 
generated. The probe was labeled with Cy5 at 5'-end and BHQ2 at 3'-end. 


Amplification 



Standard Curve 



O Standard 
X Unknown 

- Cy5 E=100.9% R A 2=1.000 Slope=3.301 y-int=38.414 


Amplification 




O Standard 
■ Unknown 

- Texas Red E=98.6% R A 2=0.995 Slope=3.357 y-int=38.119 


c The amplification curves (top) and a standard curve (bottom) for detec¬ 
tion of TGEV N gene were generated. The probe was labeled with HEX 
at 5'-end and BHQ1 at 3'-end. d The amplification curves (top) and a 
standard curve (bottom) for detection of PEDV N gene were generated. 
The probe was labeled with TexasRed at 5'-end and BHQ2 at 3'-end 
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50 nM probe for PEAV. The results demonstrated that the mul¬ 
tiplex real-time RT-qPCR could detect all target genes of these 
four viruses efficiently with high correlation values (Fig. 2). All 
the standard curves showed excellent correlation coefficient and 
amplification efficacy, for details, PEDV (R 2 = 0.999; Efif% = 
103.5), PDCoV (R 2 = 0.999; Eff%= 104.7), TGEV (R 2 = 
0.999; Eff% = 105.6), and PEAV (R 2 = 1; Eff% = 100.4) (Fig. 
2). The limit of detection of this multiplex RT-qPCR was approx¬ 
imately 100 copies of each vims per reaction (Fig. 2). 

The specificity of multiplex real-time RT-qPCR assay 

To evaluate the specificity of the multiplex RT-qPCR devel¬ 
oped in this study for swine enteric coronavims detection, the 
DNAs/cDNAs of nine other major swine vimses were used as 
templates for amplification with this multiplex system. The 
cDNAs of PEDV, PDCoV, TGEV, and PEAV were served as 
positive control while the nuclease-free water was used as a 
negative control. The results displayed that all swine enteric 
coronavimses were successfully detected. However, there was 
no positive signal detected from other nine swine vimses, 
PAstV, PKV, FMDV-O, FMDV-A, PRRSV, CSFV, PRV, 


PCV2, PrV, and negative control using this multiplex RT- 
qPCR system (Fig. 3), demonstrating that the TaqMan- 
probe-based multiplex RT-qPCR system was highly specific. 

The sensitivity of multiplex real-time RT-qPCR assay 

To determine the sensitivity of this multiplex real-time RT- 
qPCR, tenfold serial dilution of linearized plasmid mixtures 
were added to the amplification system. The results showed 
that 1.1 x 10 2 copies of TGEV, PEDV, PDCoV, and PEAV 
were detectable with the CT values at 28.75, 28.94, 30.03, 
and 28.4 respectively (Table 2). However, 1.1 x 10 1 copies 
of each target were not detectable in the same amplification 
system (Table 2). Also, the amplification exhibited reliable 
and high efficacy, PEDV (R 2 = 0.998; Eff% = 96.04), 
PDCoV (R 2 = 0.992; Eff% = 91.38), TGEV (R 2 = 0.998; 
Eff% = 94.57), and PEAV ( R 2 = 0.999; Eff% = 98.60). This 
multiplex real-time RT-qPCR demonstrated the same sensitiv¬ 
ity level when using viral RNAs as templates (Supplemental 
Fig. SI). These results implied that the TaqMan-probe-based 
multiplex real-time RT-qPCR assay could detect porcine 


Fig. 2 Establishment of multiplex 
real-time RT-qPCR. a 
Amplification curves and b stan¬ 
dard curves of optimized multi¬ 
plex TaqMan-probe-based real¬ 
time RT-qPCR for detection of 
PEDV, PDCoV, TGEV, and 
PEAV were generated at the opti¬ 
mum amplification conditions. 
The correlation coefficient and 
amplification efficacy of the stan¬ 
dards curves, PEDV ( R 2 = 0.999; 
Eff% = 103.5), PDCoV (R 2 = 
0.999; Efl% = 104.7), TGEV 
(. R 2 = 0.999; EB% = 105.6), and 
PEAV (R 2 = 1; Eff% = 100.4), 
were ideal for detecting the target 
genes 



Log Starting Quantity 

O Standard 
X Unknown 

FAM E=103.5% R A 2=0.999 Slope=3.241 y-int=35.125 

- HEX E=105.6% R A 2=0.999 Slope=3.195 y-int=34.949 

- Texas Red E=100.4% R A 2=1.000 Slope=3.312 y-int=36.136 

- Cy5 E=104.7% R A 2=0.999 Slope=3.215 y-int=35.575 
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Fig. 3 The specificity of 
multiplex TaqMan-probe-based 
real-time RT-qPCR. 1—4: positive 
templates of TGEV, PDCoV, 
PEDV, and PEAV. 5-14: tem¬ 
plates of PAsV, PKV, FMDV-O, 
FMDV-A, CSFV, PRRSV, PRV, 
PCV2, RrV, and nuclease-free 
water control 


Amplification 



0 10 20 30 40 

Cycles 


enteric coronaviruses (PEDV, TGEV, PDCoV, and PEAV) as 
less as 10 to 100 copies. 

Repeatability of multiplex real-time RT-qPCR assay 

To estimate the reproducibility of the multiplex real-time 
PCR, tenfold serial dilution of pooled linearization plasmids 
were in a triplicate manner for both intra-assay and inter-assay. 
For intra-assay, the standard plasmids were amplified three 
times simultaneously. For inter-assay, standard curves were 
done at three individual times and using a different batch of 
a standard substance. As shown in Table 3, the coefficient of 
variation for PEDV intra-assay was 0.30-1.56% and inter¬ 
assay was 0.44-1.79%. The PDCoV amplification exhibited 
the coefficients of variation at 0.09-5.51% and 1.91-3.78% 
for inter-assay and intra-assay respectively. As for TGEV, the 
coefficients of variation were 0.33-1.29% and 0.95-3.70% 
corresponding to intra-assay and inter-assay. The amplifica¬ 
tion of PEAV showed relative higher coefficients of variation, 
which was 0.39-7.32% for intra-assay and 1.08^.52% for 
inter-assay. These results indicated that the TaqMan-probe- 


Table 2 Sensitivity of multiplex RT-qPCR 


Templates (copies/gene) 

CT value 



TGEV 

PEDV 

PDCoV 

PEAV 

1.1 X 10 7 

11.54 

12.05 

12.16 

11.88 

1.1 X 10 6 

15.30 

15.66 

16.98 

15.18 

1.1 X 10 5 

18.77 

18.76 

20.37 

18.60 

1.1 X 10 4 

22.55 

22.84 

24.22 

22.38 

1.1 X 10 3 

25.71 

26.04 

27.31 

25.54 

1.1 X 10 2 

28.75 

28.94 

30.03 

28.40 

1.1 X 10 1 

ND 

ND 

ND 

ND 

NTC 

ND 

ND 

ND 

ND 


based multiplex real-time RT-qPCR assay established in this 
study was repeatable and reliable. 

Detection of the clinical samples 

The cDNAs originated from 354 clinical samples were sub¬ 
jected to amplification by this multiplex real-time RT-qPCR 
and confirmed by conventional RT-PCR (Supplemental Fig. 
S2). The results demonstrated that the samples from Henan 
and Shaanxi provinces were swine enteric coronavirus¬ 
negative while the Gansu province suffered a high prevalence 
of coronaviruses, showing 83.9% (26/31) TGEV, 48.4% (15/ 
31) PDCoV, and 38.7% (12/31) PEDV positive (Table 4). In 
Ningxia province, TGEV was also highly prevalent but the 
positive rates for the other three swine enteric coronaviruses 
were quite low or not detectable (Table 4). Only PEDV (15/ 
62, 24.2%) and PDCoV (8/62, 12.9%) among these four 
coronaviruses were detectable from clinical samples of 
Chongqing. All these four swine enteric coronaviruses were 
detectable in samples of Liaoning province, with TGEV pos¬ 
itive rate at 12.6% (18/143) and the other three were less than 
10% (Table 4). Overall, PDCoV (121/354, 34.2%) and TGEV 
(72/354, 20.3%) were the main coronaviruses detected from 
all these clinical samples. The positive rate for the newly 
emerged coronavirus in Guangdong, China, PEAV, was rarely 
found from these provinces. 

We further analyzed the co-infection of coronaviruses 
for all these positive samples. The results elucidated that 
the main co-infection were caused by a dual coronavirus, 
for example, 13 PDCoV and TGEV co-infection samples, 
five PEDV and PDCoV co-infection samples, and four 
PEDV and TGEV co-infection samples (Fig. 4). Some 
clinical samples were co-infected by three (9 PEDV, 
PDCoV, and TGEV co-infection) or all four coronaviruses 
(Fig. 4). These results indicated the pathogens of viral 
diarrhea disease in Chinese pig farms were complicated. 
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Table 3 Cq value analysis of 
multiplex real-time RT-qPCR 
repeatability 


Plasmid 

Dilution (copies/qL) 

Intra-assay 


Inter-assay 


Cq value (mean ± SD) 

CV% 

Cq value (mean ± SD) 

CV% 

PEDV 

1.1 X 10 7 

12.72 ±0.15 

1.18 

12.31 ±0.22 

1.79 


1.1 X 10 6 

16.16 ±0.05 

0.30 

15.69 ±0.07 

0.44 


1.1 X 10 5 

19.23 ±0.23 

1.20 

18.91 ±0.25 

1.32 


1.1 X 10 4 

22.85 ±0.08 

0.35 

22.67 ±0.15 

0.66 


1.1 X 10 3 

25.90 ±0.33 

1.27 

25.73 ±0.40 

1.55 


1.1 X 10 2 

28.93 ±0.45 

1.56 

28.93 ±0.37 

1.29 


1.1 X 10 1 

ND 


ND 



NTC 

ND 


ND 


PDCoV 

1.1 X 10 7 

12.69 ±0.70 

5.51 

12.71 ±0.48 

3.78 


1.1 X 10 6 

16.44 ±0.02 

0.12 

16.39 ±0.52 

3.17 


1.1 X 10 5 

19.72 ±0.60 

3.04 

19.69 ±0.60 

3.05 


1.1 X 10 4 

23.21 ±0.02 

0.09 

23.38 ±0.73 

3.12 


1.1 X 10 3 

26.14±0.52 

2.0 

26.38 ±0.82 

3.11 


1.1 X 10 2 

29.15 ±0.15 

0.51 

29.38 ±0.56 

1.91 


1.1 X 10 1 

ND 


ND 



NTC 

ND 


ND 


TGEV 

1.1 X 10 7 

11.91 ±0.07 

0.59 

11.89 ±0.44 

3.70 


1.1 X 10 6 

15.10 ±0.05 

0.33 

15.29 ±0.44 

2.88 


1.1 X 10 5 

18.32 ±0.13 

0.71 

18.74 ±0.42 

2.24 


1.1 X 10 4 

21.86 ±0.12 

0.55 

22.28 ±0.39 

1.75 


1.1 X 10 3 

24.98 ±0.32 

1.29 

25.46 ±0.26 

1.02 


1.1 X 10 2 

27.83 ±0.11 

0.40 

28.56 ±0.27 

0.95 


1.1 X 10 1 

ND 


ND 



NTC 

ND 


ND 


PEAV 

1.1 X 10 7 

12.16 ±0.89 

7.32 

12.38 ±0.56 

4.52 


1.1 X 10 6 

15.39 ±0.1 

0.65 

15.57 ± 0.51 

3.28 


1.1 X 10 5 

18.64 ±0.05 

0.27 

19.11 ±0.54 

2.82 


1.1 X 10 4 

22.14 ±0.07 

0.32 

22.60 ±0.35 

1.55 


1.1 X 10 3 

25.20 ±0.29 

1.15 

25.85 ±0.28 

1.08 


1.1 X 10 2 

28.37 ±0.11 

0.39 

28.97 ±0.57 

1.97 


1.1 X 10 1 

ND 


ND 



NTC 

ND 


ND 



Discussion 

Porcine diarrheal disease is one of the most severe dis¬ 
eases in pig farms. Swine enteric coronaviruses are the 
most significant pathogens causing diarrhea in piglets, 
especially for newborn piglets (Butler et al. 2014). All 
swine enteric coronaviruses, including traditional PEDV, 
TGEV, and newly emerged PDCoV, PEAV, could cause 
serious diarrhea, vomiting, dehydration, weight loss, and 
up to 100% death in suckling piglets (Gong et al. 2017; 
Hsu et al. 2018). Previously, some investigators devel¬ 
oped a serial of ELISA assays to detect PEDV, TGEV, 
or PDCoV infections based on M, N proteins or whole 
virus (Luo et al. 2017; Ma et al. 2016; Su et al. 2016). 


However, molecular diagnostic tools, rather than serolog¬ 
ical methods, for swine enteric coronaviruses are urgently 
needed due to their similar and high pathogenicity to 
suckling piglets, lacking a mature immune system. For 
swine enteric coronavirus detection, conventional RT- 
PCR, multiplex RT-qPCR, or pan-coronaviruses RT-PCR 
methods targeting M, N, S, or polymerase genes were 
developed in recent years (Hsu et al. 2018; Song et al. 
2006). Additionally, Zhou et al. (2018a) recently reported 
a specific real-time PCR for the detection of PEAV. To our 
knowledge, there was no such method that could simultaneous¬ 
ly detect and differentiate TGEV, PEDV, PDCoV, and PEAV 
from the same detection vial. To solve this urgent and important 
issue, we developed this multiplex real-time RT-qPCR for the 
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Table 4 Positive rates detected from the clinical samples of different provinces 
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Fig. 4 Co-infection analysis between detected swine enteric 
coronaviruses. A total of 354 clinical diarrhea samples collected from 
six provinces of China between 2015 and 2018 were subjected to an 
analysis by the real-time RT-qPCR assay developed in this study. The 
Venn diagram shows the number of samples infected by either single 
vims or multi-coronaviruses 

detection and differential diagnosis of the swine enteric 
coronaviruses circulating in pig herds. 

TaqMan-probe-based real-time qPCR is a rapid, high spec¬ 
ificity, high sensitivity, and great reproducibility detection tool 
for identification of viruses (Chang et al. 2014). Among the 
seven open reading frames (ORFs) and four structural protein 
genes, S (spike), E (envelope), M (membrane), N 
(nucleocapsid) of coronaviruses (Su et al. 2018), the M and 
N genes are highly conserved within the same antigenic group 
but less homologous between each of these four swine enteric 
coronaviruses (Yang et al. 2019). To establish a high-specific 
multiplex real-time RT-qPCR for the detection and differential 
diagnosis of swine enteric coronaviruses, primer sets, and 
TaqMan probes were designed targeting the highly conserved 
regions of PEDV or PDCoV M genes and TGEV or PEAV N 
genes, based on bioinformatics analysis of each virus. Our 
results demonstrated that each primer and probe set can only 
detect target gene itself and could not bind with any other 
targets, indicating high specificity. 

Our results showed that the singular real-time RT-qPCRs 
were capable to detect as less as 10 copies of PEDV, PDCoV, 
TGEV, and PEAV templates. However, the detection limit of 
each target gene in the multiplex real-time RT-qPCR was ap¬ 
proximately 100 copies, implying that the sensitivity of mul¬ 
tiplex real-time RT-qPCR was lower than that of singular real¬ 
time RT-qPCR probably due to the competition between 
primers, probes, templates, and reagents. 

The analysis of clinical diarrhea samples using this devel¬ 
oped multiplex RT-qPCR elucidated that the co-infection of 
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swine enteric coronaviruses commonly existed in some pig 
farms. Co-infection of swine enteric coronaviruses may cause 
recombination between co-infected viruses. During 2012 and 
2016, some new swine enteric coronaviruses were generated 
by recombination with PEDV and TGEV and spread across 
the central Eastern European countries (Akimkin et al. 2016; 
Belsham et al. 2016; Boniotti et al. 2016). The recombination 
might create more virulent enteric virus strains or new viruses, 
leading to potential outbreaks or pandemics of swine viral 
diarrhea. Additionally, co-infections may promote the evolu¬ 
tion of non-pathogenic enteric viruses into high virulent and 
pathogenic viruses. A recent example showed the existence of 
PDCoV in pig herds for many years in China mainland and 
Hong Kong but without observed clinical symptoms (Pan 
et al. 2017; Woo et al. 2012). Unfortunately, this virus caused 
severe outbreaks in some states of the USA since 2014 (Ma 
et al. 2015; Wang et al. 2014). This virulence change was most 
probably caused by the previous exposure to PEDV or other 
co-infected swine enteric viruses. Therefore, swine enteric 
coronaviruses co-infection may bring risks for the prevention 
and control of swine diarrhea. 

In summary, we developed a TaqMan-probe-based multiplex 
real-time RT-qPCR with high specificity and sensitivity for si¬ 
multaneous detection and differential diagnosis of swine enteric 
coronaviruses, PEDV, TGEV, PDCoV, and PEAV. This real-time 
RT-qPCR assay is of great significance for the prevention and 
control and epidemiological investigation of swine viral diarrhea. 
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